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We derive a continuum equation for the magnetization of a conducting ferromagnet in the presence 
of a spin-polarized current. Current effects enter in the form of a topological term in the Landau- 
Lifshitz equation . In the stationary situation the problem maps onto the motion of a classical 
charged particle in the field of a magnetic monopole. The spatial dependence of the magnetization 
is calculated for a one-dimensional geometry and suggestions for experimental observation are made. 
We also consider time-dependent solutions and predict a spin-wave instability for large currents. 
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Phenomena associated with spin-polarized currents in 
layered materials and in Mn-oxides have attracted high 
interest recently. Efforts are strongly concentrated on 
theoretical and experimental investigation of large mag- 
netoresistance, which is of great value for future applica- 
tions. Examples of the effect are GMR in layered materi- 
als (see review [0), spin valve effect for a particular case 
of a three-layer sandwich and CMR in the manganese 
oxides (see review B). 

The dependence of resistivity on magnetic field is ex- 
plained conceptually in two steps: first the magnetic field 
changes the magnetic configuration of the material and 
that in turn influences the current. Of course, as for any 
interaction there must be a back-action of the current 
on the magnetic structure. The existence of such back- 
action was explored in |||]]. Several current-controlled 
micro-devices utilizing this principle were proposed 
In both papers layered structures with magnetization be- 
ing constant throughout the magnetic layers were consid- 
ered. In the present paper we derive the equations for a 
continuously changing magnetization in the presence of 
a spin-polarized current. This equation takes the form of 
a Landau-Lifshitz equation with an additional topologi- 
cal term, and admits a useful analogy with a mechanical 
system. We discuss several solutions in one-dimensional 
geometries. Our equations also can be viewed as a con- 
tinuum generalization of ||,[|] for layer thickness going to 
zero. 

Consider a current propagating through a conducting 
ferromagnet. Conducting electrons are viewed as free 
electrons interacting only with local magnetization M. 
The motion of each individual electron is governed by 
the Schroedinger equation with a term Jf/erM, where 
Jh is the value of the Hund's rule coupling or in gen- 
eral of the local exchange. Since spin-up electrons have 
lower energy a nonzero average spin of conducting elec- 
trons (1/2) (er) develops. An angular momentum density 
(h/2)(er) is then carried with the electron current so we 
have a flux of angular momentum. This leads to a non- 
zero average torque acting on the magnetization which 



can deflect it from the original direction (see figure 1). 
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FIG. 1. Experimental setting: spin-polarized current en- 
ters a half-infinite magnet from the left. Originally the mag- 
netization is aligned along the easy-axis v\\z. However if the 
incoming electrons are spin-polarized in a different direction, 
their interaction with the magnetization leads to a deflection 
of the magnetization. 

Propagation of a current in a ferromagnet should be de- 
scribed by a system of two equations: one for the motion 
of conducting electrons and another for the magnetiza- 
tion. We derive here the second equation in the limit of 
small space-time gradients and present several solutions. 
The case of very large Jh — > oo is considered, meaning a 
complete polarization of electron spins in the ferromag- 
net. This case can be often realized in experiment. In 
the layered structures magnetic layers can be made of a 
material with large band splitting, like Heusler alloys, in 
which the spin opposite to magnetization direction can 
not propagate. In the CMR materials large Hund's rule 
coupling is well known Q and constitutes the basis for 
a double-exchange mechanism governing their magnetic 
ordering. 

Schroedinger Equation: the conducting electrons 
are considered noninteracting with an e = p 2 /2m energy 
spectrum: 
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We diagonalize the matrix M(r, t)cr a p with a local spin 
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rotation 4> a = U a p( r,t)<ij)p. The spinor <j> describes the 
electron in the coordinate system with z-axis being par- 
allel to the local magnetization. Retaining only the first 
order terms in gradients and using 0_ — > for Jh — > oo 
we reduce (flf) from a system of two equations to one 



equation for 



ih-zr = A(b 

dt 2m 



spin amplitude: 
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The last term in (0) can be transformed (M = Mn 
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where A mon (n) is a function satisfying the following 
equations: 
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If we view n 2 (r,t) = 1 as a sphere, A mon has the sim- 
ple interpretation of the vector potential due to a mag- 
netic monopole located at the center of the sphere. The 
monopole term is known to appear from U + pV iUp + in 
the theory of Berry phase and is used by other CMR 
theories in different forms (see ||). Equation (||) has a 
form of a Schroedinger equation in a magnetic field ex- 
panded up to the linear term in A off , with vector poten- 
tial A° s = (ih 2 /4m)A^ non \7 i n k . It describes the motion 
of the conducting electrons in the given field n(r, t). Con- 
versely it gives the interaction between the current and 
the magnetization. The form of equation is the same 
as for an electromagnetic interaction, and hence we can 
write by analogy: 
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where j is an electric current. 

Magnetization motion is described by Landau- 
Lifshitz equations which are obtained from the energy 
functional. After adding (|]) to the usual energy density 
of a ferromagnet with uniaxial anisotropy along the axis 
v , we obtain: 
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with K > corresponding to easy-axis and K < to 
easy-plane magnets. The equations of motion then take 
the form: 
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where the last term in f is new and describes the effect 
of the current. The system of (||) and (^,0) constitute 
a complete set of equations for a magnet with current. 
Equations (^|J^), generalizing the Landau-Lifshitz equa- 
tion in the presence of a current, are the central result of 
this work. 

Since magnetization corresponds to angular momen- 
tum L = H/gfiB M, an equation of the angular momen- 
tum flux continuity follows from (^): 
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The flux Afci consists of two parts: one due to the spatial 
derivatives of magnetization and another due to the mo- 
tion of conducting electrons. In our situation the spins 
of moving electrons are parallel to n. That is why their 
contribution is factorized in the form (h/2)(ji/e)rik- 

Consider the stationary case in an experimental 
setting shown on figure 1 . For the stationary process the 
r.h.s. of equation ([|) vanishes. The current propagates 
along the y direction. All spatial derivatives reduce to 
V — > Vj,. For the reasons immediately following we will 
denote differentiation with a prime to get a resemblance 
to a time derivative in notation V y n = n'. From (Q) we 
get an equation on n(r, t): 

J[n" x n] = [(-K(vn)u - Q(-)W x nl) x n] (10) 

e 

with new parameters: 
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Since n in the stationary case depends on y only, we can 
interpret y a s a fictitious time; together with n 2 = 1 



equation (10) can then be interpreted as the equation of 
motion for a particle of a mass J confined to the surface 
of a unit sphere and experiencing two forces: 

(a) a force of magnitude — K(yw)v parallel to the 
anisotropy axis, 

(b) a Lorentz force , due to a field H mon = — Q(^)n of a 
magnetic monopole positioned in the center of the sphere. 
The vector product ensures that only tangential compo- 
nents of the total force act on the particle. The normal 
component is compensated by the reaction forces. Such 
an anology enables one to visualize the soluctions of the 
original equation ( |l0| ) as trajectories of a massive particle 
on the sphere. 

The equation of particle motion in the field of a mag- 
netic monopole ([h]) has two first integrals M. 
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Together they give a way to solve (|l0| ) for arbitrary ini- 
tial conditions. Expressing everything through the Euler 
angles {(f>(y), 9(y)} (defined on figure 1) of the vector n , 
we obtain: 

, D v - Q(j/e) cos 9 



J sin 2 6 



Ql _ 2W-K cos 2 9 (D v - Q(j/e) cos 9) 2 

~ V J J 2 sin 2 9 

= F{9) (13) 

The problem for 9 is solved by the implicit function: 
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afterwards <f> can be found from the first equation in (13). 

Assume that deep inside the magnet (y — > oo) the 
magnetization resumes its original direction along the 
anisotropy axis n — > v, n' — > 0. From this the val- 
ues of the first integrals can be found and substituted 
into (|l3|). Natural length and current scales appear in 
the calculation: 
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through which the material parameters J, K, M enter the 
problem. Their values for different materials are given in 
the following table: 



Material 


L m [A] 


jo [A/cm 2 ] 


Ref. 


CMR: Lao.eeCao.asMnOa 


> 130 


< 4 • 10 y 






Fe 


40 


1.1 • 10 8 






Heusler Alloy: PtMnSb 


50 - 100 


- 5 • 10 Y 





The integral (|1J) can be then expressed in elementary 
functions but the formula is long and will be detailed in 
a later paper. Instead, the results are presented on fig- 
ure 2. It is seen that magnetization relaxes in a distance 
~ 10L m which is about the width of the domain wall in 
the material. 

In the "particle picture" the motion starts at some 
point A on the trajectory, yet to be determined from the 
b.c. on the normal metal - magnet interface, and ends 
on the North pole. The particle has just enough energy 
to climb the potential hill and come to rest on the top. 
The particle trajectory is bent by the monopole field. In 
the absence of the monopole the particle would go along 
the meridian. 

The boundary condition on the metal-magnet 
interface, y = is derived from the continuity of the 
angular momentum flux. Such condition ensures that 
there is no torque concentrated on the boundary consis- 
tent with the assumption of slow spatial changes of the 
magnetization. 




FIG. 2. Magnetization deflection angle 8(y). Curves corre- 
spond to different values of current, 6o = t/2, p = 0.5. Inset: 
trajectory of particle on the sphere for j — jo- Starting point 
A corresponds to the conditions of the main graph. 

The reflection of the down-spin electron component oc- 
curs on the length scale of the electron wavelength. On 
this distance magnetization is almost constant and solv- 
ing the one-particle reflection problem we find the jump 
of the electron flux component in the y direction Ej = £ 
to be: 
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where £ is the average injected flux. From (||) the flux 
inside the magnet is: 



A yiclcctron 

l ^i-y ~ ■"magnet 



M 2 J[n x n' ]i 
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In the metal, only the electron part of the flux is present. 
Then continuity gives the boundary condition: 



S - (Sn )n = M 2 J[n x n' ] 



(18) 



Note that it involves both the vector n and its derivative 
on the boundary. 

Condition ([18]) can be transformed into a system of 
two algebraic equations and an inequality: 



p(cos#o — xy) = 1 — x 

y 2 + x 2 — 2yx cos 9q < sin 2 9q 



(19) 



where S = £ e, cos 9q = e • v, x = v ■ no, y = e • no- 
The parameter p = eYl/j, p E [0,1/2], describes the 
"degree of polarization" of the incident electrons . The 
inequality in ( |l9| ) is a geometrical constraint on x and y 
arising from their definition. 

The trajectory is determined by three parameters: 
(j/jo,p,9o). We can plot a domain of existence of a 
solution to (jl^) in the 3-D space of these parameters. 
A typical 2-D section of this diagram for constant 9q is 
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shown in figure 3. A solution is absent in the regions 
B and C which means that for larger currents and spin- 
polarizations no smooth stationary solution approaching 
the easy-axis direction at infinity is available. Either 
a non-stationary solution or a solution which never ap- 
proaches v will be realized in that region. 
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FIG. 3. A typical 2-D section of the phase diagram plotted 
for 6o = it/3. A: domain of existence of a solution to (|l9|); B: 
no solution, b.c. at y = can not be satisfied; C: no solution, 
b.c. at y — > oo can not be satisfied and spin wave instability 
occurs. 

Time-dependent solutions of equation (0) can be 
found in some cases. We again assume the current j to 
be uniform. We rewrite equation (0) through n(r,t): 



^ = [g x n] - Q(^)V in 
at e 

g = JAn + K(un)v 



(20) 



and suppose no(r) solves [g x n ] = 0, i.e. represents a 
static solution in the absence of the current. Then 
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n(r,t) = n (r + Q(-)t) = n (r + ^-J , 
e V2 Jo 



t) (21) 



where ujq — K, is a solution of (|20|) for a nonzero current. 
For instance a moving Bloch wall will be a solution when 
current is flowing perpendicular to it (provided pinning 
is absent). 

Another particular solution is a spin wave in the pres- 
ence of a current. We search for a solution (Q) in the 
form of a spin wave: {8 — const, <j) = kr — ujt}. This 
gives the spectrum: 



uo = Jk 2 -Q 



K 



(22) 



As we see, the current changes the energy gap of spin 
waves and shifts the position of the minimum: 
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where a is the angle between j and k and luq = K is 
the gap of spin wave in an anisotropic ferromagnet. For 
large enough current j > jo an instability occurs. 
That is also the condition which leads in the region C 
on figure 3 to the loss of any trajectory approaching v 



at infinity as the integral (|14j) becomes undetermined. A 
spin-wave instability is also predicted in other models of 
spin-polarized transport p| . 

Discussing possible experiments we note that the 
characteristic current is large, but such densities are in 
fact common for layered metallic structures and j ~ jo 
is experimentally possible. In this regime the calcu- 
lated magnetization profile (figure 2) shows a deviation of 
« 20° on the boundary. Detection of the effect is difficult 
because the spatial resolution of existing probes exceeds 
L m . For a quantitative measurement of the deflection 
angle on the boundary, element specific X-ray magnetic 
circular dichroism (MXCD) JllJ could be used. A sin- 
gle layer of different magnetic element can be put on the 
boundary in the process of the film growth. Such a layer 
will not much disturb the overall magnetization profile. 
The MXCD signal of the additional layer could be sepa- 
rated from the Mn signal and thus (9(0) can be measured. 
The time-dependent approach developed here can be ap- 
plied to other experimental geometries, e.g. to find the 
continuum analogies for the devices proposed in || . 
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